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Abstract 

We consider the fifth order partial differential equation (PDE) 

U4:x,t - 5uxxt + ^ut + uu^x + 2uxUix - huu^x - lQuxU.j,x + I2uu.j, = 0, 

which is a generalization of the integrable Camassa-Holm equation. 
The fifth order PDE has exact solutions in terms of an arbitrary num- 
ber of superposed pulsons, with geodesic Hamiltonian dynamics that 
is known to be integrable in the two-body case N ~ 2. Numerical sim- 
ulations show that the pulsons are stable, dominate the initial value 
problem and scatter elastically. These characteristics are reminiscent of 
solitons in integrable systems. However, after demonstrating the non- 
existence of a suitable Lagrangian or bi-Hamiltonian structure, and 
obtaining negative results from Painleve analysis and the Wahlquist- 
Estabrook method, we assert that the fifth order PDE is not integrable. 

1 Introduction 

This note is concerned with the fifth order partial differential equation 
(PDE) 

UAx,t - ^Uxxt + ^Ut + UUbx + '^UxU4x - ^UU^x - lOUxUxx + 12uUx = 0. (1.1) 

One reason for our interest in this equation is that it admits exact solutions 
of the form 

N 

u = ^pj-(t)(2e-l^-5^WI - e-2|^-*(*)l), (1.2) 
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where pj , qj satisfy the canonical Hamiltonian dynamics generated by 

1 ^ 

Hn = -Y. PjPfc(2e-'*-«*' - e-2|*-5^l). (1.3) 
j,k=i 

Following lH], we refer to such solutions as "pulsons." The equations for 
the A^-body pulson dynamics are equivalent to geodesic flow on an A^- 
dimensional space with coordinates qj and co-metric 

g'' = 9iQj-qk), <?(x)=2e-l-'l-e-2N. 

The pulsons ( |l.2D are weak solutions with discontinuous second derivatives 
at isolated points. 
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Figure 1: Pulson solutions (|1.2| ) of equation ( |l.lD emerge from a Gaussian 
of unit area and width a = 5 centered about x = 33 on a periodic domain 
of length L = 100. The fastest pulson crosses the domain four times and 
collides elastically with the slower ones. 



The PDE (LI) is one of a family of integral partial differential equations 
considered in given by 

rrit + urux + 2uxm = 0, u = g * m, (1-4) 

where u{x, t) is defined in terms of m{x, t) by the convolution integral 

j-OO 

g*rn:= g{x - y)m{y,t) dy. 
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The integral kernel g{x) is taken to be an even function, and for any g the 
equation ( |1.4| ) has the Lie-Poisson Hamiltonian form 

6H 

mt = -{mdx + d^m)— (1.5) 

where 

If If 
H= — / mg*mdx = — / mudx. (1-6) 

Any equation in this family admits pulson solutions 

N 

u{x,t) = ^pj{t)g{x - qj{t)) 

for arbitrary N, with Pj,qj satisfying the canonical Hamilton's equations 

dpj dHN // N dqj SHn ^ , . 

-dt = = "^^ f^^'^" " ^ = ^ = f^^'^^^^^- " 

(1.7) 

generated by the Hamiltonian 
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The equations ( |1.7| ) correspond to geodesic motion on a manifold with co- 
metric g^^ = g{qj — q^)- A significant result of ||l| is that the two-body 
dynamics {N = 2) is integrable for any choice of kernel g, and numerical 
calculations show that this elastic two-pulson scattering dominates the initial 
value problem. 

Three special cases are isolated in Q, namely (up to suitable scaling) 

• g{x) = 6{x) - Riemann shocks, 

• g{x) = 1 — l^l, |x| < 1 - compactons, 

• g{x) = e~l^l - peakons. 



For each of these cases both the integral PDE ( |1.4| ) and the corresponding 
finite-dimensional system ( |1.7D (for any A'^) are integrable. Of most relevance 
here is the third case, where g{x) = e"'^', which is the (scaled) Green's 
function for the Helmholtz operator, satisfying the identity 

il-dl)g{x) = 26ix). 

In that case after rescaling we may take 

m = u-Uxx, (1-8) 
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Figure 2: Two rear-end collisions of pulson solutions ( |1.2| ) of equation ( |1.1| ). 
The initial positions are x = 25 and x = 75. The faster pulson moves at 
twice the speed of the slower one. For this ratio of speeds, both collisions 
result in a phase shift to the right for the faster space-time trajectory, but 
no phase shift for the slower one. 



and the equation ( [I.4| ) is just a PDF for u{x,t), namely 

Ut - Uxxt - UU3x - 2UxUxx + SuUx = 0, (1.9) 

which is the dispersionless form of the integrable Camassa-Holm equation 
for shallow water waves |Q, ^] . For Camassa-Holm the pulson solutions take 
the form of peakons or peaked solitons, i.e. 

N 

=^Pj(t)e-l^-*WI. (1.10) 
i=i 

The fifth order equation arises from a different choice of Green's function. 
Using the identity 

(4 - dl){l - dl)g{x) = m{x), g{x) = 2e-l-l - e'^l-l (1.11) 

we find that (after suitable scaling) this choice of g yields 

m = — 5uxx + 4ti, (1-12) 



and then the equation (1^) becomes the fifth order equation (IT). Thus 
the PDF ( |1.1| ) should be regarded as a natural higher order generalization 
of the Camassa-Holm equation. 

2 Hamiltonian and Lagrangian considerations 

In a forthcoming article [Q] we discuss a more general family of integral 
PDEs of the form 

mt + urux + buxTn = 0, u = g * m, (2T) 
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where b is an arbitrary parameter; ( |1.4| ) corresponds to the particular case 
6 = 2, and numerical results have recently been established for different b 
values in B. In the case of the peakon kernel g = e"'^', with m given by 



(1.8), the equations in this class were tested by the method of asymptotic 
integrability and only the cases 6 = 2,3 were isolated as potentially 
integrable. For 6 = 2 the integrability of the Camassa-Holm equation by 
inverse scattering was already known |2|, ^, but for the new equation 

Ut - Uxxt - UU3x - "iUxUxx + ^UU^ = 0, (2.2) 

with 6 = 3 the integrability was proved in 0] by the construction of the 
Lax pair. The two integrable cases 6 = 2,3 were also found recently via the 
perturbative symmetry approach ||8|. For any 6 7^ —1, the peakon family 
(2.1) with g = e~l^l arises as the dispersionless limit at the quadratic order 



in the asymptotic expansion for shallow water waves Q. 

Another motivation for our interest in the fifth order equation (^]^) is 



that it is expressed naturally in terms of the quantity m (1.12), given by the 
product of two Helmholtz operators acting on u. Such a product appears in 
the fifth order operator 

Bo = dx{A-dl){l-dl), 

which we found Q to provide the first Hamiltonian structure for the the new 
integrable equation (2.2). This led us to the conjecture that the operator Bq 



should appear naturally in the theory of higher order integrable equations 
such as ( ll.ll) . 

All the equations of the form ( |l.4D have the Lie-Poisson Hamiltonian 
structure given by (|1.5| ), but for integrability we expect a bi-Hamiltonian 
structure. In the case of the Camassa-Holm equation ( |l.g| ) there are two 
ways to derive a second Hamiltonian structure. The first is by inspection 
using a conservation law, noting that (|1.9p may be written as 



mt = yuuxx + -^u^ - J = dx-^ = dx{l - d^)— (2.3) 

for 

H = -\ j{uul + u^)dx. (2.4) 

The identity ( p.3| ) gives the second Hamiltonian structure for Camassa- 
Holm, and mdx + dxm, 3^(1 — 3^) constitute a compatible bi-Hamiltonian 
pair. 

Similarly j mdx is conserved for (1.1), with m given by (1.12). The 
conservation law is explicitly 

{uix-^Uxx+^u)t = - (^uu4x + UxUsx " ^^L " ^'^i'^^xx " + =: ^x- 

' (2.5) 
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By analogy with the Camassa-Holm equation, this would suggest that a 
suitable constant coefficient Hamiltonian operator might be 5a;(4— 3^) 
(which we know to be a Hamiltonian operator for the new equation ( |2.2[ )). 
This would require the right hand side of ( |2.5| ) to take the form 



However, for the flux of (|2.5|) we find 



ou 

for any local density functional K of u, and we suppose that the operators 
dx{4: — — d'^) and mdx + dxm must be incompatible. 

The second way to derive the Hamiltonian structure ( |2.3| ) for Camassa- 
Holm is via the action (integral of Lagrangian density) 

S = j j C[4>]dxdt:= j j ^ {(t>x(t>t - Mt + (l)x4>L + </'x) dx dt, 
for u = (px- A Legendre transformation yields the conjugate momentum 

dH 1 TTl 

^=2('A.-'/'3.) = y, 
and the same Hamiltonian as (12.41) above, i.e. 



^ = j Q"^</'t - dx. 
Trying the same approach for ( |1.1D , we set u = (j)x and rewrite it as 

(t>hx,t - 54>3x,t + Hxt + (t)x4'Gx + 2(^a:x<A5x " 5<px4>4x " W(t)xx4>3x + 'i-^(t>x4>xx = 0. 

(2.6) 

However, the equation (|2.6| ) cannot be derived from a local Lagrangian den- 
sity due to the presence of the terms (px^'Gx + '^4'xx4'5x- 

The first nonlocal Hamiltonian structure for the Camassa-Holm equation 
is obtained by applying the recursion operator to mdx + dxm. This means 
that ( ll.Sp can be written in the Hamiltonian form 

mt = {mdx + dxm){dx - dx)~^{mdx + ^^m) — , H = mdx. 

With the same H, the analogous identity for (LI) is 

mt = B— = {mdx + dxm){dl - f>dl + 45,)"^ (m9, + dxm) — , (2.7) 
dm dm 

but from the above considerations we would expect that the formal nonlocal 
operator B on the right hand side of ( |2.7D is not Hamiltonian, and indeed 
using the functional equations derived in Q it is possible to show that it 
fails to satisfy the Jacobi identity. 
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3 Reciprocal transformation and Painleve analysis 



Having failed to find the sort of Lagrangian or bi-Hamiltonian structure for 



(1.1) that we would reasonably expect, we proceed to see what Painleve 
analysis can tell us about this fifth order equation. However, we note that 
both the Camassa-Holm equation ( |1.9| ) and the new equation ( |2.2| ) provide 
examples of the weak Painleve property |jl^ , with algebraic branching in the 
solutions. For these equations we have found it convenient to use reciprocal 



transformations (see |11] for definitions), which transform to equations with 
pole singularities, and indeed in Q this was the key to our discovery of the 
Lax pair for (|2.2|). Hodograph transformations of this kind have been used 



before to remove branching from classes of evolution equations |12, |l^], but 
here we are dealing with non-evolutionary PDEs. 

To make the results of our analysis more general, we will consider the 
whole class of equations 

rrit + urrix + buxtn = 0, m = u/^x — 5uxx + 4u, b ^ 0, (3.1) 



for arbitrary nonzero b, which is the particular family of equations (2.1 



corresponding to the integral kernel ( 1.11 ), and includes (LI) in the special 



case 6 = 2. Each equation in the class ( |3.1| ) has the conservation law 

{m'/% = -{m^"u)x, 

and so introducing a new dependent variable p according to 

= -m (3.2) 

means that we may consistently define a reciprocal transformation to new 
independent variables X, T given by 

dX = pdx — pudt, dT = dt. (3-3) 

Transforming the derivatives we have the new conservation law 

{p-^)T = ux. (3.4) 



Rewriting the relation ( 1.12 ) in terms of dx and using ( p.4| ) to eliminate 



derivatives of u we obtain the identity 

u = \{f>- {pdx?) ■ ipdx) ■P{p-')T - J, (3.5) 
which means that ( |3.4| ) can be written as an equation for p alone, i.e. 

(p'')t = (J {(pdx? - 5)p{logp)xT - . (3.6) 
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The fifth order equation ( |3.6| ) is the reciprocal transform of ( p.ip . Rather 
than carrying out the fuh Painleve test for the transformed equation, it is 
sufficient for our purposes to follow [14| and apply the Painleve test for 



ODEs to the travelling wave reduction of ( |3.6|) . Hence we set p = p{z), 
z = X — cT and the resulting fifth order ODE may be integrated twice to 
get the third order ODE 

sUJ "n^'^ -^^2-^)-W' = w^)^~P^"' 

(3.7) 

6 7^ 1, where d, e are arbitrary constants and c is the wave speed, with prime 
denoting d/dz. For 6=1 there is a logp term on the right hand side, and 
so this case has logarithmic branching and is immediately excluded by the 
Painleve test. Similarly, because of the p^~^ term all non-integer values of h 
have branching and are discarded. 

We proceed to apply Painleve analysis to (|3.7|) for integer h ^ 0, 1, 
seeking leading order behaviour at a movable point zq of the form p ~ 
a{z — zq)^ for integer exponent /i. For all integers h < —2 the only possible 
balance is fi = 4/(3 — 6), which is non- integer and hence gives algebraic 
branching. In the special case 6 = — 1 there are four possible balances with 
= 1, with a and the resonances depending on the value of c; we have 
checked that no value of c gives all integer resonances, so the Painleve test 
is failed. For the remaining cases of integer 6 > 2 we find n = 1 with a? = \ 
or = 4. For integer 6 > 4 there is also the balance ^ = 4/(3—6) which is in 
general non-integer, except for the special cases /i = —4, —2, —1 for 6 = 4, 5, 7 
respectively. Thus all integer values of 6 apart from 6 = 2,3, 4, 5, 7 are ruled 
out by the (strong) Painleve test due to algebraic branching; but they could 
still be analysed by the weak Painleve test if we allow such branching. 

Let us consider in more detail the first two types of balance for integer 
6 > 2. When p ~ ib(z — zq) we have a non-principal balance with resonances 
r = -1,-1,3. Interestingly, the resonance condition at r = 3 is failed 
when 6 = 2, the obstruction being the term (so for these balances 
the test is only passed in the limit c ^ oo), but satisfied for all integer 
6 > 3. However, for the principal balances p ~ ib2(z — zq) the resonances 
are r = —1,1/2,3/2 which means there is algebraic branching and so the 
(strong) Painleve test is failed for any 6. We have further checked whether 



the weak Painleve test of [10| could be satisfied by allowing an expansion 
in powers of {z — zqY^'^ in the principal balance. However, the resonance 
condition is satisfied at r = 1/2 but failed at r = 3/2, meaning that this 
expansion with square root branching cannot represent the general solution 
as it doesn't contain enough arbitrary constants. The arbitrariness can 
only be restored by adding infinitely many terms in powers of log(2; — zq), 
and so no form of Painleve property can be recovered. The existence of 
logarithmic branching in both the principal and non-principal balances is a 
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strong indication of non-integrability. 

It is interesting to observe that when the first term on the right hand 
side of ( |3.7| ) is absent (the hmit c — > oo), it admits exact solutions in terms 
of trigonometric/hyperbolic functions, corresponding to the first order re- 
ductions 

{p'f = l + 2dp + \{^e- (f)p^, (p'f = A + 8dp + ^{2d^ - e)/. 

In fact we can also see that the original equation (|l.l| ) fails the weak 
Painleve test directly. For the Camassa-Holm equation ( |1.9| ) the test is 
satisfied by a principal balance 

2 

u ~ -(l)t/4>x + a(p3 + ... 
with resonances —1,0,2/3, with the singular manifold 4>{x,t) and a{x,t) 



being arbitrary. For (1.1) there is an analogous balance 
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u ~ -(pt/<Px + a(p3 + 

with resonances -1,0,4/3, (l±V41)/6; the presence of irrational resonances 
implies logarithmic branching. 



4 Prolongation algebra method 

While Painleve analysis is a good heuristic tool for isolating potentially 
integrable equations, it can never be said to provide definite proof of non- 
integrability. If one gives a precise definition of integrability in terms of 
existence of infinitely many commuting symmetries, then the symmetry 
approach of Shabat et al [^] gives necessary conditions for integrability 
(but does not provide a constructive way to find a Lax pair or linearization 
when such conditions are satisfied). The symmetry approach has only very 
recently been extended Q so that it can be applied to nonlocal or non- 
evolution equations such as (pTll), (ll.9|) . As an alternative, we apply the 
prolongation algebra method of Wahlquist and Estabrook |16|, and directly 



seek a Lax pair for (1.1) in the form of a compatible linear system 

^x = U^, % = (4.1) 

for suitable matrices U, V (usually taking values in the fundamental repre- 
sentation of a semi-simple Lie algebra) which should depend on u and its 
derivatives, and on a spectral parameter. We have found the clear presen- 
tation of the method in [^] very useful. 



The compatibility of the system (4.1) yields the zero curvature equation 



Ut - Vx + [[/, V] = 0, (4.2) 
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and the essence of the Wahlquist-Estabrook method is that given the orig- 
inal PDE (in this case (|1.1|)) one may use (4^) to derive the functional 
dependence of U, V on ti, u^-, etc. A negative result means that no Lax pair 
of a suitable form exists, suggesting that the equation is not integrable, but 
of course this is sensitive to the initial assumptions that are made on the 
functional form of U^V . 

For ease of notation we will denote the nth derivative Unx = ^n- Given 
that ( |1.1D can be written as a conservation law for m as in (|2.5|), a reasonable 
ansatz is to assume that 

U = U{m), V = V{u,ui,U2,U3,Ui) 

(with dependence on the spectral parameter suppressed). Given the known 
form of the zero curvature representations for the equations (p.. 9]), ( |2.2D we 
further assume that U is linear in m, so that 

U = Am + B = {u4- 5u2 + ^u)A + B, 

where A, B are constant matrices (independent of x, i, but potentially de- 
pendent on the spectral parameter). Substituting this ansatz into (^^), and 
using ( ^.5| ) to eliminate the t derivative m^, we find 

(-litis - 2uiU4 -I- Sutis -I- 10uiU2 - 12^1)^ - Us - '"4K3 - 'W3K2 - '"2K1 

-uiVu + (U4 - 5u2 + 4u) [A, V] + [B,V]={) (4.3) 

(with subscripts on V denoting partial derivatives). None of the matrices 
depend on M5, so ([4.3| ) is linear in u^. In particular, the coefficient of 
must vanish, giving the equation Vu^ = —uA which integrates immediately 
to yield 

V = -uuiA + T{u,ui,U2,u^), (4.4) 

where F is so far arbitrary and must be determined from the remaining 
terms in (4.3). 

At the next step we substitute for V in the rest of ( [4.3| ) to obtain 

(— ■U1U4 -I- 5UU3 -I- lQuiU2 — \2uUi)A — U^Tu-i — ""31^112 ~ U2Tui — UiTu 

+ (ii4 - 5ii2 -h4'u)[A,F] +uui[A,B] + [B,T] = 0. (4.5) 
The coefficient of M4 gives the equation 

F„3 = [A, F + uB] - uiA 

which can be integrated exactly as 

F = e''3^A(u, ui, U2)e"''^^ - uiu^^A - uB, (4.6) 
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where A is the arbitrary function of integration. From (^]^) we see the 
presence of Adexpu^A = exp^adu^A) acting on A, which would imply 
exponential-type dependence on in the Lax pair unless {adusA)^A = 
for some positive integer n. Such exponential dependence would seem 
unlikely given that the original equation (1.1) is polynomial in u and its 
derivatives, and we will seek assumptions that prohibit infinitely many non- 
zero commutators occurring in ( [4.6[ ). 

Substituting for F from ( ^ ) in the tt4-independent terms of ( ^ ) and 
applying Ad exp(— us^) we get 

{u2U3 + 5uU3 + 10uiU2-12uUi)A-U3Au2 -U2A„^ -UiAu + {-5U2+4:U)[A, A] 

+ [e-"^^5e"^^, A] + (mug + 5uu2 - V)e-"3^Ce"3^ = 0, (4.7) 
where we have set 

C = [A,B]. 



Potentially (4.7) is an infinite power series in 113, each coefficient of which 
must vanish. The simplest assumption we can make to terminate the series 
is to take 

[A,C]=0 (4.8) 

which implies 

Ade-^'^^iB) = e"'^(-"3A)(B) = B-U3C, ^de""«^(C) = e"'='(-"3^) (C) = C, 

and hence ( [4.7D becomes linear in U3. A fortunate consequence of (^^) is 
that the coefficient of gives 

Ana = (^2 + 5u)A - [C, A] + UlC, 

which integrates exactly without further assumption to yield 

A = e-"2^-E(u, ui)e"2^ + Qn^ + 5uu2^ A + U1U2C. (4.9) 
The remaining terms in (|4.7| ), after acting with Ad exp U2C, now become 
{5uiU2-12uui)A- (^^ul + 4u'^^ C - U2Eu^ - uiEu + {-^U2 + 4.u)[A, E\ 

Once again we are faced with an infinite power series, this time in U2- Before 
looking for further simplifying assumptions, we note that the coefficient of 
the term linear in U2 is just 

-Eu, + ui{bA + [B, C]) - [5A + [B, CIE]= 0, 
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which integrates immediately to 



D = bA+[B,C]. 



(4.11) 



To analyse the other terms in ( 4.10| ) we find it convenient to introduce the 
quantities 

F=[C,[C,B\], G=[D,Bl 



and note that the identities 
[A,D]=Q = [A,F], [A,G] 



[C,D]=F, [[B,C],D]=f) (4.12) 



aU hold. 

The coefficient of u'2 in (P^ ) is then 

-lc+\[F,E]-uiF = Q, 
and the coefficient of U2 gives 

-12nni^ - 4n2e"i^Ce~"i^ - mZ^ 



(4.13) 



+ ^u[A, Z] + 



Be 



Z + -uiD 
2 ^ 



(4.14) 



(after substituting for E from ( [l.llD and acting with Ade^^^). We shall not 
need to consider the equations [{adCy^B,E] — nui{adC)^B = occurring 
at U2 , re > 3. Instead we look at the coefficient of ui in ( 4.14 ), which is 

-Zu + [G, Z] - Au^F - UuA = 0. (4.15) 

We are unable to integrate this directly without making a further assump- 
tion, the simplest possible being 



which implies 



[F, G] = 



3 



(4.16) 



(4.17) 



We must now use the remaining equations to determine the commutation 
relations for the constant Lie algebra elements A, B, C, D, F, G, Q. 
Looking at the coefficient of Ui in ( [4.13 ) we see that 

-^C+^-[F,Z]=0, 

and using (|1|) and ( p^ with (|1|) yields 

[F,e] = 3e~^"C7e^", 
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which immediately imphes 



[F,e]=3C7, [C,G]=0. (4.18) 

Using the Jacobi identity we also have 

= [C, G] = [C, [D, B]] = -[B, [C, D]] - [D, [B, C]] = [B, F], 

which by further applications of the Jacobi identity gives 

[C,F]=0 = [D,F]. (4.19) 

Now we return to the equation ( [4.13 ) and use ( [4.11 ) to evaluate the coeffi- 
cient of ui as 

-^[F,[D,Z]]-F = 0. (4.20) 

Substituting for Z as in ( [4.17[ ) and taking the constant coefficient in 
([420| ) we use ([4l2| ), (|1|), (|1D to find 

= -^[F, [D, e]]-F= i ([9, [F, D]] + [D, [6, F]]) - F 



--[D,C]-F 



'-F 
2 



Then F = implies C = from the first equation in (4.18), and it is 
straightforward to show that the Lax pair (4.1) collapses down to the trivial 
case [U, V] = 0, with ([4.2|) reducing to the scalar equation nit = as in 



5 Conclusions 

Since the fifth order PDE ([l.l[ ) is in the class of pulson equations studied 
by Pringer and Holm it admits exact solutions in the form of a direct 
superposition of an arbitrary number of pulsons (as in Figure |l]). These 
particular A^-pulson solutions take the precise form ([1.2[) , and by the general 
results of |l|] we know that for N = 2 the Hamiltonian equations ( [1.7[ ) 
describing the two-body dynamics are integrable. However, several different 



considerations provide strong evidence that the fifth order PDE (1.1) is not 
integrable in the sense of admitting a Lax pair and being solvable by the 
inverse scattering transform. 

It is well known that integrable PDEs such as the Korteweg-deVries 
equation or Camassa-Holm ||2[ ^ admit a compatible pair of Hamiltonian 
structures which together define a recursion operator generating infinitely 
many higher symmetries. We have tried unsuccessfully to find an analogous 
bi-Hamiltonian or Lagrangian formulation for the fifth order equation ([1.1[), 
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but as far as we are aware it admits only the single Hamiltonian structure 
d). 

Both the Camassa-Holm equation ( |1.9| ) and the new integrable equa- 
tion ( ^.21) isolated by Degasperis and Procesi Q exhibit the weak Painleve 
property of , with algebraic branching in local expansions representing a 
general solution. There are many examples of Liouville integrable systems 
in finite dimensions |18] and Lax integrable PDEs with this property. 
For evolution equations transformations of hodograph type can restore the 
strong Painleve property |l^, and similarly in [Q, ^ we have used recipro- 
cal transformations for the non-evolutionary equations (^]^) , ( |2.2| ) . In order 
to apply Painleve analysis to (^]^) we have found it convenient to employ 
a reciprocal transformation which removes the branching at leading order, 
but further analysis of the travelling wave reduction shows that there is still 
algebraic branching in the principal balances due to half-integer resonances. 
Furthermore, in both principal and non-principal balances a resonance con- 
dition is failed, and so after the transformation even the weak Painleve test 
cannot be satisfied. 

We have also applied an integrability test which is perhaps less fash- 
ionable nowadays, namely the prolongation algebra method of Wahlquist 
and Estabrook ||l^ . By making certain simple assumptions we find that no 
polynomial Lax pair of a suitable form exists for ( |1.1| ). 

We cannot expect the iV-body pulson dynamical system to be integrable 
for arbitrary N > 2, since this would imply the existence of an infinite 
integrable subsector within a non-integrable PDF. However, it would be 
good to find an analytical explanation for the apparent soliton-like behaviour 
of the pulson solutions, and their numerical stability as evidenced by Figure 
^. Further analytical and numerical studies will be required to understand 
the stability properties of the pulson solutions. 
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